THE KINETIC THEORY OF DILUTE SOLUTIONS OF FLEXIBLE 
POLYMERS: HYDRODYNAMIC INTERACTION 

J. Ravi Prakash a 

a Department of Chemical Engineering, Indian Institute of Technology, 
Madras, India, 600 036 

1. INTRODUCTION 

The rheological properties of dilute polymer solutions are commonly used in industry for characterising 
the dissolved polymer in terms of its molecular weight, its mean molecular size, its chain architecture, the 
relaxation time spectrum, the translational diffusion coefficient and so on. There is therefore considerable 
effort world wide on developing molecular theories that relate the microscopic structure of the polymer and 
its interactions with the solvent to the observed macroscopic behavior. In this chapter, recent theoretical 
progress that has been made in the development of a coherent conceptual framework for modelling the 
rheological properties of dilute polymer solutions is reviewed. 

A polymer solute molecule dissolved in a dilute Newtonian solvent is typically represented in molecular 
theories by a coarse-grained mechanical model, while the relatively rapidly varying motions of solvent 
molecules surrounding the polymer molecule are replaced by a random force field acting on the mechanical 
model. The replacement of the complex polymer molecule with a coarse-grained mechanical model is 
justified by the belief that such models capture those large scale properties of the polymer molecule, such 
as its stretching and orientation by the solvent flow field, that are considered to be responsible for the 
solution's macroscopic behavior. An example of a coarse-grained model frequently used to represent a 
flexible polymer molecule is the bead-spring chain, which is a linear chain of identical beads connected 
by elastic springs. 

Progress in the development of molecular theories for dilute polymer solutions has essentially involved 
the succesive introduction, at the molecular level, of various physical phenomena that are considered to be 
responsible for the macroscopic properties of the polymer solution. For instance, the simplest theory based 
on a bead-spring model assumes that the solvent influences the motion of the beads by exerting a drag 
force and a Brownian force. Since this theory fails to predict a large number of the observed features of 
polymer solutions, more advanced theories have been developed which incorporate additional microscopic 
phenomena. Thus, theories have been developed which (i) include the phenomenon of 'hydrodynamic 
interaction' between the beads, (ii) try to account for the finite extensibility of the polymer molecule, 
(iii) attempt to ensure that two parts of the polymer chain do not occupy the same place at the same 
time, (iv) consider the internal friction experienced when two parts of a polymer chain close to each 
other in space move apart, and so on. The aim of this chapter is to present the unified framework within 
which these microscopic phenomena may be treated, and to focus in particular on recent advances in the 
treatment of the effect of hydrodynamic interaction. To a large extent, the notation that is used here is 
the same as that in the treatise Dynamics of Polymeric Liquids by Bird and co-authors [ ^| . 

2. TRANSPORT PROPERTIES OF DILUTE SOLUTIONS 

2.1. Dilute solutions 

A solution is considered dilute if the polymer chains are isolated from each other and have negligible 
interactions with each other. In this regime of concentration the polymer solution's properties are deter- 
mined by the nature of the interaction between the segments of a single polymer chain with each other, 
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and by the nature of the interaction between the segments and the surrounding solvent molecules. As the 
concentration of polymers is increased, a new threshold is reached where the polymer molecules begin to 
interpenetrate and interact with each other. This threshold is reached at a surprisingly low concentration, 
and heralds the inception of the semi-dilute regime, where the polymer solution's properties have been 
found to be significantly different. Beyond the semi-dilute regime lie concentrated solutions and melts. 
In this chapter we are concerned exclusively with the behavior of dilute solutions. 

A discussion of the threshold concentration at which the semi-dilute regime in initiated is helpful in 
introducing several concepts that are used frequently in the description of polymer solutions. 

A polymer molecule surrounded by solvent molecules undergoes thermal motion. A measure of the 
average size of the polymer molecule is the root mean square distance between the two ends of the polymer 
chain, typically denoted by R. This size is routinely measured with the help of scattering experiments, 
and is found to increase with the molecular weight of the polymer chain with a scaling law, R ~ M v , 
where M is the molecular weight, and v is the scaling exponent which depends on the nature of the 
polymer-solvent interaction. In good solvents, solute-solvent interactions are favoured relative to solute- 
solute interactions. As a consequence the polymer chain swells and its size is found to scale with an 
exponent v = 3/5. On the other hand, in poor solvents, the situation is one in which solute-solute 
interactions are preferred. There exists a particular temperature, called the theta temperature, at which 
the scaling exponent v changes dramatically from 3/5 to 1/2. At this temperature, the urge to expand 
caused by two parts of the chain being unable to occupy the same location (leading to the presence of an 
excluded volume), is just balanced by the repulsion of the solute molecules by the solvent molecules. 

Polymer chains in a solution can be imagined to begin to interact with each other when the solution 
volume is filled with closely packed spheres representing the average size of the molecule. This occurs 
when ri v R 3 w 1, where n p in the number of chains per unit volume. Since n p = p p N\/M, where p p is the 
polymer mass density and Na is Avagadro's number, it follows that polymer density at overlap, p*, scales 
with molecular weight as, p* ~ M Y ~ Zv . Polymer molecules typically have molecular weights between 10 4 
and 10 6 gm/mol. As a result, it is clear that the polymer solution can be considered dilute only at very 
low polymer densities. Since experimental measurements are difficult at such low concentrations, the 
usual practice is to extrapolate results of experiments carried out at decreasing concentrations to the 
limit of zero concentration. For instance, in the case of dilute polymer solutions it is conventional to 
report the intrinsic viscosity, which is defined by, 

[rj\ = lim -HL- 
p p ^o p p r\ s 

where r\ p is the polymer contribution to the solution viscosity, and r\ s is the solvent viscosity. 

2.2. Homogeneous flows 

Complex flow situations typically encountered in polymer processing frequently involve a combination 
of shearing and extensional deformations. The response of the polymer solution to these two modes of 
deformation is very different. Consequently, predicting the rheological properties of the solution under 
both shear and extensional deformation is considered to be very important in order to properly charac- 
terise the solutions behavior. Rather than considering flows where both these modes of deformation are 
simultaneously present, it is common in polymer kinetic theory to analyse simpler flow situations called 
homogeneous flows, where they may be treated separately. 

A flow is called homogeneous, if the rate of strain tensor, 7 = (Vv)(t) + (Vv)t(i), where v is the solution 
velocity field, is independent of position. In other words, the solution velocity field v in homogeneous 
flows, can always be represented as v = v + n(t) ■ r, where v is a constant vector, /e(i) = Vv(t) is a 
traceless tensor for incompressible fluids, and r is the position vector with respect to a laboratory fixed 
frame of reference. While there is no spatial variation in the rate of strain tensor in homogeneous flows, 
there is no restriction with regard to its variation in time. Therefore, the response of dilute solutions 
to transient shear and extensional flows is also used to probe its character as an alternative means of 
characterisation independent of the steady state material functions. 
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Two homogeneous flows, steady simple shear flow and small amplitude oscillatory shear flow, that 
are frequently used to validate the predictions of molecular theories which incorporate hydrodynamic 
interaction, are described briefly below. A comprehensive discussion of material functions in various flow 
situations can be found in the book by Bird et at [ . 

2.3. Simple shear flows 

The rheological properties of a dilute polymer solution can be obtained once the stress tensor, r, is 
known. The stress tensor is considered to be given by the sum of two contributions, r = t s + t p , where 
t s is the contribution from the solvent, and r p is the polymer contribution. Since the solvent is assumed 
to be Newtonian, the solvent stress (using a compressive definition for the stress tensor [ |]) is given by, 
t s = — r) s 7. The nature of the polymer contribution t p in simple shear flows is discussed below. 

Simple shear flows are described by a velocity field, 

Vx = lyx V, Vy = 0,V Z =0 (2) 

where the velocity gradient ^ yx can be a function of time. From considerations of symmetry, one can 
show that the most general form that the polymer contribution to the stress tensor can have in simple 
shear flows is [ , 

\ 

(3) 

where the matrix of components in a Cartesian coordinate system is displayed. The form of the stress 
tensor implies that only three independent combinations can be measured for an incompressible fluid. 
All simple shear flows are consequently characterised by three material functions. 

2.3.1. Steady simple shear flows 

Steady simple shear flows are described by a constant shear rate, 7 = ^yx | . The tensor k is consequently 
given by the following matrix representation in the laboratory-fixed coordinate system, 

(4) 

The three independent material functions used to characterize such flows are the viscosity, ri P , and the 
first and second normal stress difference coefficients, ^1 and 1 i>2, respectively. These functions are defined 
by the following relations, 

T%y = "7 Vp ; T P X - T P y = - 7 2 *! ; T P y - T P Z = -j 2 * 2 (5) 

where r p y , T p x ,r p y are the components of the polymer contribution to the stress tensor t p . 

At low shear rates, the viscosity and the first normal stress coefficient are observed to have constant 
values, rjpfi and 'i'l,!), termed the zero shear rate viscosity and the zero shear rate first normal stress 
coefficient, respectively. At these shear rates the fluid is consequently Newtonian in its behavior. 

At higher shear rates, most dilute polymer solutions show shear thinning behavior. The viscosity 
and the first normal stress coefficient decrease with increasing shear rate, and exhibit a pronounced 
power law region. At very high shear rates, the viscosity has been observed to level off and approach a 
constant value, r) p ooi called the infinite shear rate visosity. A high shear rate limiting value has not been 
observed for the first normal stress coefficient. The second normal stress coefficient is much smaller in 
magnitude than the first normal stress coefficient, however its sign has not been conclusively established 
experimentally. Note that the normal stress differences are zero for a Newtonian fluid. The existence of 
non-zero normal stress differences is an indication that the fluid is viscoelastic. 
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Experiments with very high molecular weight systems seem to suggest that polymer solutions can also 
shear thicken. It has been observed that the viscosity passes through a minimum with increasing shear 
rate, and then increases until a plateau region before shear thinning again [ [T^], 

It is appropriate here to note that shear flow material functions are usually displayed in terms of 
the reduced variables, rj p /rj p .Q, ^E'i/ v I / i l o and ^2/^1, versus a non-dimensional shear rate /?, which is 
defined by f3 = A P 7, where, A p = [r)]o M rj s /Np^ k& T, is a characteristic relaxation time. The subscript 
on the square bracket indicates that this quantity is evaluated in the limit of vanishing shear rate, 
fee is Boltzmann's constant and T is the absolute temperature. For dilute solutions one can show that, 
[v]/[v]o = Vp/Vpfi and (3 = <rj P;0 i/n p k B T. 

2.3.2. Small amplitude oscillatory shear flow 

A transient experiment that is used very often to characterise polymer solutions is small amplitude 
oscillatory shear flow. The upper plate in a simple shear experiment is made to undergo sinusoidal 
oscillations in the plane of flow with frequency 10. For oscillatory flow between narrow slits, the shear 
rate at any position in the fluid is given by [ |J, jyxit) — 70 cosojt, where 70 is the amplitude. The tensor 
n(t) is consequently given by 

/0 1 0\ 

n(t) = 70 cos ut (6) 
\0 0/ 

Since the polymer contribution to the shear stress in oscillatory shear flow, t^ x , undergoes a phase shift 
with respect to the shear strain and the strain rate, it is customary to represent its dependence on time 
through the relation [ |l[ , 

t£ x = —rj'(uj) 70 cos cut — r)"(uj) j sin uit (7) 

where rf and rf' are the material functions characterising oscillatory shear flow. It is common to represent 
them in a combined form as the complex viscosity, rf = rf — i rf' . 

Two material functions which are entirely equivalent to rf and 77" and which are often used to display 
experimental data, are the storage modulus G' = cur]" / '(nksT) and the loss modulus G" = wf|'/(nfeT). 
Note that the term involving G' in equation(^) is in phase with the strain while that involving G" is in 
phase with the strain rate. For an elastic material, G" — 0, while for a Newtonian fluid, G' = 0. Thus, 
G 1 and G" are measures of the extent of the fluid's viscoelasticity. 

In flow situations which have a small displacement gradient, termed the linear viscoelastic flow regime, 
the stress tensor in polymeric fluids is described by the linear constitutive relation, 

dsG(t - s)j(t,s) (8) 

> 

where G(t) is the relaxation modulus. 

When the amplitude 70 is very small, oscillatory shear flow is a linear viscoelastic flow and consequently 
can also be described in terms of a relaxation modulus G(t) . Indeed, expressions for the real and imaginary 
parts of the complex viscosity can be found from the expression, 

y* 00 

rf=\ G(s)e' tujs ds (9) 
Jo 

Experimental plots of log G' and log G" versus nondimensional frequency show three distinct power 
law regimes. The regime of interest is the intermediate regime [ where for dilute solutions of high 
molecular weight polymers in good or theta solvents, both G' and G" have been observed to scale with 
frequency as a; 2 / 3 . 
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It is appropriate to note here that the zero shear rate viscosity 77 Pj o & n d the zero shear rate first normal 
stress difference "J 1,0, which arc linear viscoclastic properties, can be obtained from the complex viscosity 
in the limit of vanishing frequency, 



r) p>0 = lim n'(w) ; 



*i,o 



lim 

u—>0 



(10) 



2.4. Scaling with molecular weight 

We have already discussed the scaling of the root mean square end-to-end distance of a polymer 
molecule with its molecular weight. In this section we discuss the scaling of the zero shear rate intrinsic 
viscosity [rj\o, and the translational diffusion coefficient D, with the molecular weight, M. As we shall 
see later, these have proven to be vitally important as experimental benchmarks in attempts to improve 
predictions of molecular theories. 

It has been found that the relationship between [77] and M can be expressed by the formula, 



where, a is called the Mark-Houwink exponent, and the prefactor K depends on the polymer-solvent 
system. The value of the parameter a lies between 0.5 and 0.8, with the lower limit corresponding to 
theta conditions, and the upper limit to a good solvent with a very high molecular weight polymer solute. 
Measured intrinsic viscosities are routinely used to determine the molecular weight of samples once the 
constants K and a are known for a particular polymer-solvent pair. 

The translational diffusion coefficient D for a flexible polymer in a dilute solution can be measured by 
dynamic light scattering methods, and is found to scale with molecular weight as [ |2j, 



where the exponent fj, lies in the range 0.49 to 0.6. Most theta solutions have values of /i close to the lower 
limit. On the other hand, there is wide variety in the value of /i reported for good solvents. It appears 
that the upper limit is attained only for very large molecular weight polymers and the intermediate 
values, corresponding to a cross over region, are more typical of real polymers with moderate molecular 
weights. 

2.5. Universal behavior 

It is appropriate at this point to discuss the most important aspect of the behavior of polymer solu- 
tions (as far as the theoretical modelling of these solutions is concerned) that is revealed by the various 
experimental observations. When the experimental data for high molecular weight systems is plotted in 
terms of appropriately normalized coordinates, the most noticeable feature is the exhibition of universal 
behavior. By this it is meant that curves for different values of a parameter, such as the molecular weight, 
the temperature, or even for different types of monomers can be superposed onto a single curve. For 
example, when the reduced intrinsic viscosity, [t?]/[t?]o is plotted as a function of the reduced shear rate 
P, the curves for polystyrene in different types of good solvents at various temperatures collapse onto a 
single curve [ [l| . 

There is, however, an important point that must be noted. While polymers dissolved in both theta 
solvents and good solvents show universal behavior, the universal behavior is different in the two cases. 
An example of this is the observed scaling behavior of various quantities with molecular weight. The 
scaling is universal within the context of a particular type of solvent. The term universality class is used 
to describe the set of systems that exhibit common universal behavior [ |l| . Thus theta and good solvents 
belong to different universality classes. 

The existence of universality classes is very significant for the theoretical description of polymer so- 
lutions. Any attempt made at modelling a polymer solution's properties might expect that a proper 
description must incorporate the chemical structure of the polymer into the model, since this determines 
its microscopic behavior. Thus a detailed consideration of bonds, sidegroups, etc. may be envisaged. 



Mo = KM' 



a 



(ii) 



D ~ M'f 



(12) 



5 



However, the universal behavior that is revealed by experiments suggests that macroscopic properties of 
the polymer solution are determined by a few large scale properties of the polymer molecule. Structural 
details may be ignored since at length scales in the order of nanometers, different polymer molecules 
become equivalent to each other, and behave in the same manner. As a result, polymer solutions that 
differ from each other with regard to the chemical structure or molecular weight of the polymer molecules 
that are dissolved in it, the temperature, and so on, still behave similarly as long as a few parameters 
that describe molecular features are the same. 

This universal behavior justifies the introduction of crude mechanical models, such as the bead-spring 
chain, to represent real polymer molecules. On the other hand, it is interesting to note that in many 
cases, the predictions of these models are not universal, ft turns out that apart from a basic length and 
time scale, there occur other parameters that need to be prescribed, for example, the number of beads N 
in the chain, the strength of hydrodynamic interaction h*, the finite spring extensibility parameter b, and 
so on. ft is perhaps not incorrect to state that any molecular theory that is developed must ultimately 
verify that universal predictions of transport properties are indeed obtained. The universal predictions 
of kinetic theory models with hydrodynamic interaction are discussed later on in this chapter. 



3. BEAD-SPRING CHAIN MODELS 

The development of a kinetic theory for dilute solutions has been approached in two different ways. 
One of them is an intuitive approach in the configuration space of a single molecule, with a particular 
mechanical model chosen to represent the macromolecule, such as a freely rotating bead-rod chain or a 



freely jointed bead-spring chain [ |13|, [56|, 45 1. The other approach is to develop a formal theory in the 
phase space of the entire solution, with the polymer molecule represented by a general mechanical model 
that may have internal constraints, such as constant bond lengths and angles [ [ll| H ||. The results of 
the former method are completely contained within the latter method, and several ad hoc assumptions 
made in the intuitive treatment are clarified and placed in proper context by the development of the 
rigorous phase space theory. Kinetic theories developed for flexible macromolecules in dilute solutions 
have generally pursued the intuitive approach, with the bead-spring model proving to be the most popular. 
This is because the lack of internal constaints in the model makes the formulation of the theory simpler. 
Recently, Curtiss and Bird [ acknowledging the 'notational and mathematical' complexity of the 
rigorous phase space theory for general mechanical models, have summarised the results of phase space 
theory for the special case of bead-spring models with arbitrary connectivity, ie. for linear chains, rings, 
stars, combs and branched chains. 

In this section, since we are primarily concerned with reviewing recent developments in theories for 
flexible macromolecules, we describe the development of kinetic theories in the configuration space of a 
single molecule. However, readers who wish the understand the origin of the ad hoc expressions used 
for the Brownian forces and the hydrodynamic force, and the formal development of expressions for the 
momentum and mass flux, are urged to read the article by Curtiss and Bird [ || . 

The general diffusion equation that governs the time evolution of the distribution of configurations 
of a bead-spring chain subject to various nonlinear effects, and the microscopic origin of the polymer 
contribution to the stress tensor are discussed in this section. The simplest bead-spring chain model, the 
Rouse model is also discussed. We begin, however, by describing the equilibrium statistical mechanical 
arguments that justify the representation of a polymer molecule with a bead-spring chain model, and we 
discuss the equilibrium configurations of such a model. 



3.1. Equilibrium configurations 

When a flexible polymer chain in a quiescent dilute solution is considered at a lowered resolution, ie. 
at a coarse-grained level, it would appear like a strand of highly coiled spaghetti, and the extent of its 
coiling would depend on its degree of flexibility. A quantity used to characterise a chain's flexibility is the 
orientational correlation function, whose value K or (At) , is a measure of the correlation in the direction 
of the chain at two different points on the chain which are separated by a distance At along the length 



G 



of the chain. At sufficiently large distances A£, it is expected that the correlations vanish. However, it is 
possible to define a persistence length £ ps , such that for A£ > ^ ps , orientational correlations are negligible 

[©■ 

The existence of a persistence length suggests that as far as the global properties of a flexible polymer 
chain are concerned, such as the distribution function for the end-to-end distance of the chain, the 
continuous chain could be replaced by a freely jointed chain made up of rigid links connected together at 
joints that are completely flexible, whose linear segments are each longer than the persistence length £ ps , 
and whose contour length is the same as that of the continuous chain. 

The freely jointed chain undergoing thermal motion is clearly analogous to a random-walk in space, 
with each random step in the walk representing a link in the chain assuming a random orientation. Thus 
all the statistical properties of a random-walk are, by analogy, also the statistical properties of the freely 
jointed chain. The equivalence of a polymer chain with a random-walk lies at the heart of a number of 
fundamental results in polymer physics. 

3.1.1. Distribution functions and averages 

In polymer kinetic theory, the freely jointed chain is assumed to have beads at the junction points 
betwen the links, and is referred to as the freely jointed bead-rod chain [ 0]. The introduction of the 
beads is to account for the mass of the polymer molecule and the viscous drag experienced by the polymer 
molecule. While in reality the mass and drag are distributed continuously along the length of the chain, 
the model assumes that the total mass and drag may be distributed over a finite number of discrete 
beads. 

For a general chain model consisting of N beads, which have position vectors r V) v = 1, 2, . . . , N, in a 
laboratory fixed coordinate system, the Hamiltonian is given by, 

H = tC + ^(r 1 ,r 2 ,...,r N ) (13) 

where K, is the kinetic energy of the system and <j) is the potential energy. 4> depends on the location of 
all the particles. 

The center of mass r c of the chain, and its velocity r c are given by 
1 N 1 N 

v=l v=l 

where r u = dr v jdt. The location of a bead with respect to the center of mass is specified by the vector 

If Qi, Q2, ■ ■ ■ Qd denote the generalised internal coordinates required to specify the configuration of 
the chain, then the kinetic energy of the chain in terms of the velocity of the center of mass and the 
generalised velocities Q s = dQ s /dt, is given by [ |^], 

K = ^^ + i]T5> t Q s Q, (15) 

s t 

where the indices s and t vary from 1 to d, m is the mass of a bead, and g s t is the metric matrix, defined 
by, g s t — m J2 U {dR v / dQ s ) ■ (dR v / dQt)- In terms of the momentum of the center of mass, p c = mNr c , 
and the generalised momenta P s , defined by, P s — (dlC/dQ s ), the kinetic energy has the form [ 0, 

IC = — !— p 2 r + - V V G st Ps Pt (16) 

s t 

where, G s t are the components of the matrix inverse to the metric matrix, ^\ G s t gtu = <^s«, and 5 SU is 
the Kronecker delta. 
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The probability, V cq dr c dQ dp c dP, that an TV-bead chain model has a configuration in the range 
dr c dQ about r c , Q and momenta in the range dp c dP about p c , P is given by, 

V cq (r c , Q,p c , P) =Z- 1 e- n / k * T (17) 
where Z is the partition function, defined by, 

Z = 1 1 1 1 e- n ^ kBT dr c dQdp c dP (18) 



The abbreviations, Q and dQ have been used to denote Qi, Q2, ■ ■ ■ , Qd and dQi dQ2 ■ ■ ■ dQd , respec- 
tively, and a similar notation has been used for the momenta. 

The configurational distribution function for a general iV-bead chain, ip cq (Q ) dQ, which gives the 
probability that the internal configuration is in the range dQ about Q , is obtained by integrating V cq 
over all the momenta and over the coordinates of the center of mass, 

ip cq (Q)=Z- 1 JJJ e- n ' k * T dr c dp c dP (19) 

For an TV-bead chain whose potential energy does not depend on the location of the center of mass, the 
following result is obtained by carrying out the integrations over p c and P [ 0] , 



MQ) = J7M^^dQ (20) 

where, g(Q) = det(g st ) = 1/ det(C? s4 ). 

An expression that is slightly different from the random-walk distribution is obtained on evaluating 
the right hand side of equation ( p0| ) for a freely jointed bead-rod chain. Note that the random-walk 
distribution is obtained by assuming that each link in the chain is oriented independently of all the other 
links, and that all orientations of the link are equally likely. On the other hand, equation (|2(]) suggests 
that the probability for the links in a freely jointed chain being perpendicular to each other, for a given 
solid angle, is slightly larger than the probability of being in the same direction. Inspite of this result, 
the configurational distribution function for a freely jointed bead-rod chain is almost always assumed 
to be given by the random- walk distribution [ 0]. Here afterwards in this chapter, we shall refer to a 
freely jointed bead-rod chain whose configurational distribution function is assumed to be given by the 
random- walk distribution, as an ideal chain. For future reference, note that the random-walk distribution 
is given by, 




ip eq (9i, . . . ,0 N ~i,4>i, . . . ,0 N -i) = \ — \ sin (9i (21) 



where &i and <j>i are the polar angles for the ith link in the chain [ |) . 

Since the polymer chain explores many states in the duration of an observation quantities observed on 
macroscopic length and time scales are averages of functions of the configurations and momenta of the 
polymer chain. A few definitions of averages are now introduced that are used frequently subsequently 
in the chapter. 

The average value of a function X ( r c , Q, p c , P ) , defined in the phase space of a polymer molecule is 
given by, 

(X) eq = J J J J XV cq dr c dQd Pc dP (22) 

We often encounter quantities X that depend only on the internal configurations of the polymer chain and 
not on the center of mass coordinates or momenta. In addition, if the potential energy of the chain does 
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not depend on the location of the center of mass, then it is straight forward to see that the equilibrium 
average of X is given by, 



(X> eq = J X^dQ (23) 

3.1.2. The end-to-end vector 

The end-to-end vector r of a general bead-rod chain can be found by summing the vectors that represent 
each link in the chain, 

N-l 

r = V" aui (24) 

i=l 

where a is the length of a rod, and Ui is a unit vector in the direction of the ith link of the chain. 
Note that the components of the unit vectors Ui, i = 1,2,.. . ,N — 1, can be expressed in terms of the 
generalised coordinates Q [ || . 

The probability -P e q(f) dr, that the end-to-end vector of a general bead-rod chain is in the range dr 
about r can be found by suitably contracting the configurational distribution function ipeq (Q) [ ||] , 



p cq(r) = J S^r - ^ flM,j ?/> oq ( Q 



) dQ (25) 



where 5(.) represents a Dirac delta function. 

With V'cq ( Q ) given by the random- walk distribution (|2l]) , it can be shown that for large values of N 
and r = \r\ < 0.5Na, the probability distribution for the end-to-end vector is a Gaussian distribution, 




P ^ r)= { 2«{N-lW ) eXP [ 2(N-l)a> ) (26) 
The distribution function for the end-to-end vector of an ideal chain with a large number of beads N is 



therefore given by the Gaussian distribution (26) 



The mean square end-to-end distance, ( r 2 ) e , for an ideal chain can then be shown to be, ( r 2 ^ = 

(N — l)a 2 . This is the well known result that the root mean square of the end-to-end distance of a 
random-walk increases as the square root of the number of steps. In the context of the polymer chain, 
since the number of beads in the chain is directly proportional to the molecular weight, this result implies 
that R ~ M - 5 . We have seen earlier that this is exactly the scaling observed in theta solvents. Thus one 
can conclude that a polymer chain in a theta solvent behaves like an ideal chain. 

3.1.3. The bead-spring chain 

Consider an isothermal system consisting of a bead-rod chain with a constant end-to-end vector r, 
suspended in a bath of solvent molecules at temperature T. The partition function of such a constrained 
system can be found by contracting the partition function in the constraint-free case, 



For an TV-bead chain whose potential energy does not depend on the location of the center of mass, the 
integrations over r c , p c and P can be carried out to give, 



Z(r) = C J 5\r-Y^au^jit> eq (Q)dr c dQ 



(28) 
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Comparing this equation with the equation for the end-to-end vector fl25[), one can conclude that, 

Z(r) =CP eq (r) (29) 

In other words, the partition function of a general bead-rod chain (except for a multiplicative factor inde- 
pendent of r) is given by P eq (r). This result is essential to establish the motivation for the introduction 
of the bead-spring chain model. 

At constant temperature, the change in free energy accompanying a change in the end-to-end vector 
r of a bead-rod chain, by an infinitesimal amount dr, is equal to the work done in the process, ie., 
dA = F ■ dr, where F is the force required for the extension. The Helmholtz free energy of a general 
bead-rod chain with fixed end-to-end vector r can be found from equation (|29|), 

A(r) = - k B T lnZ(r) = A - k B T lnP cq (r) (30) 

where A Q is a constant independent of r. For an ideal chain, it follows from equations ( pgf ) and (|3(]), that 
a change in the end-to-end vector by dr, leads to a change in the free energy dA, given by, 

dA ^=w?h r - dr ^ 

Equation ( |3l| ) implies that there is a tension F in the ideal chain, F — (3k B T '/ (N — l)a 2 ) r, which resists 
any attempt at chain extension. Furthermore, this tension is proportional to the end-to-end vector r. 
This implies that the ideal chain acts like a Hookean spring, with a spring constant H given by, 

*-W%p < 32 > 

The equivalence of the behavior of an ideal chain to that of a Hookean spring is responsible for the 
introduction of the bead-spring chain model. Since long enough sub-chains within the ideal chain also 
have normally distributed end-to-end vectors, the entire ideal chain may be replaced by beads connected 
to each other by springs. Note that each bead in a bead-spring chain represents the mass of a sub-chain 
of the ideal chain, while the spring imitates the behavior of the end-to-end vector of the sub-chain. 

The change in the Helmholtz free energy of an ideal chain due to a change in the end-to-end vector is 
purely due to entropic considerations. The internal energy, which has only the kinetic energy contribution, 
does not depend on the end-to-end vector. Increasing the end-to-end vector of the chain decreases the 
number of allowed configurations, and this change is resisted by the chain. The entropic origin of the 
resistance is responsible for the use of the phrase entropic spring to describe the springs of the bead-spring 
chain model. 

The potential energy S, of a bead-spring chain due to the presence of Hookean springs is the sum of 
the potential energies of all the springs in the chain. For a bead-spring chain with N beads, this is given 

by 

N-l 

S = -Hj2Qi'Qi ( 33 ) 
»=i 

where Q i = ri + \ — r,; is the bead connector vector between the beads i and i + 1. The configurational 
distribution function for a Hookean bead-spring chain may be found from equation ( pp| ) by substituting 
4>(Q) — S, with the Cartesian components of the connector vectors chosen as the generalised coordinates 
Q s . The number of generalised coordinates is consequently, d = 3N — 3, reflecting the lack of any 
constraints in the model. Since g(Q) is a constant independent of Q for the bead-spring chain model [ 
H, one can show that, 




^ oq (Qi,...,Q J v-i)= 7T-r^ exph— -Q 7 .0- (34) 
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It is clear from equation ( p4| ) that the equilibrium distribution function for each connector vector in 
the bead-spring chain is a Gaussian distribution, and these distributions are independent of each other. 
From the property of Gaussian distributions, it follows that the vector connecting any two beads in a 
bead-spring chain at equilibrium also obeys a Gaussian distribution. 

The Hookean bead-spring chain model has the unrealistic feature that the magnitude of the end-to- 
end vector has no upper bound and can infact extend to infinity. On the other hand, the real polymer 
molecule has a finite fully extended length. This deficiency of the bead-spring chain model is not serious 
at equilibrium, but becomes important in strong flows where the polymer molecule is highly extended. 
Improved models seek to correct this deficiency by modifying the force law between the beads of the 
chain such that the chain stiffens as its extension increases. An example of such a nonlinear spring force 
law that is very commonly used in polymer literature is the finitely extensible nonlinear elastic (FENE) 
force law [ || . 

3.1.4- Excluded volume 

The universal behavior of polymers dissolved in theta solvents can be explained by recognising that all 
high molecular weight polymers dissolved in theta solvents behave like ideal chains. However, a polymer 
chain cannot be identical to an ideal chain since unlike the ideal chain, two parts of a polymer chain 
cannot occupy the same location at the same time. In the very special case of a theta solvent, the 
excluded volume force is just balanced by the repulsion of the solvent molcules. In the more commonly 
occuring case of good solvents, the excluded volume interaction acts between any two parts of the chain 
that are close to each other in space, irrespective of their distance from each other along the chain length, 
and leads to a swelling of the chain. This is a long range interaction, and as a result, it seriously alters the 
macroscopic properties of the chain. Indeed there is a qualitative difference, and this difference cannot 
be treated as a small perturbation from the behavior of an ideal chain [ ^| . Curiously enough however, 
all swollen chains behave similarly to each other, and modelling this universal behavior was historically 
one of the challenges of polymer physics [ ^j, @, ||, [|. Here, we very briefly mention the manner in 
which the problem is formulated in the case of bead-spring chains. 

The presence of excluded volume causes the polymer chain to swell. However, the swelling ceases when 
the entropic retractive force balances the excluded volume force. The retractive force arises due to the 
decreasing number of conformational states available to the polymer chain due to chain expansion. This 
picture of the microscopic phenomenon is captured by writing the potential energy of the bead-spring 
chain as a sum of the spring potential energy and the potential energy due to excluded volume interactions. 
The excluded volume potential energy is found by summing the interaction energy over all pairs of beads 
/x and u, E — (1/2) ^^=1 E (r v — r„) , where E (r v — tv) is a short-range function usually taken as, 

E (r„ — r M ) = vkftT5{r v — r M ) ; v being the excluded volume parameter with dimensions of volume. 
The total potential energy of a Hookean bead-spring chain with (5-function excluded volume interactions 
is consequently, 

JV-l N 

( t >= 2 H £ Qi-Qi+2 VknT £ Hrv-r^) (35) 

i = l M.^=l 

The equilibrium configurational distribution function of a polymer chain in the presence of Hookean 



springs and excluded volume can be found by substituting equation ( p5| ) into equation (20), and all 
average properties of the chain can be found by using equation (^). Solutions to these equations in the 
limit of long chains have been found by using a number of approximate schemes since an exact treatment 
is impossible. The most accurate scheme involves the use of field theoretic and renormalisation group 
methods [ ^| . The universal scaling of a number of equilibrium properties of dilute polymer solutions with 
good solvents are correctly predicted by this theory. For instance, the end-to-end distance is predicted 
to scale with molecular weight as, R ~ M - 588 . 

The spring potential in equation (^5|) has been derived by considering the Helmholtz free energy of an 
ideal chain, ie. under theta conditions. It seems reasonable to expect that a more accurate derivation of 
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the retractive force in the chain due to entropic considerations would require the treatment of a polymer 
chain in a good solvent. This would lead to a non-Hookean force law between the beads [[?], ^9|. Such non- 
Hookean force laws have so far not been treated in non-equilibrium theories for dilute polymer solutions 
with good solvents. 

3.2. Non-equilibrium configurations 

Unlike in the case of equilibrium solutions it is not possible to derive the phase space distribution 
function for non-equilibrium solutions from very general arguments. As we shall see here it is only possible 
to derive a partial differential equation that governs the evolution of the configurational distribution 
function by considering the conservation of probability in phase space, and the equation of motion for 
the particular model chosen. The arguments relevent to a bead-spring chain are developed below. 

3.2.1. Distribution functions and averages 

The phase space of a bead-spring chain with N beads can be chosen to be given by the 6N — 6 
components of the bead position coordinates, and the bead velocities such that, 



V (ri, .. . ,r N , fx,. .. ,r N , t) dri . .. dr N dr\... dr 



N 



is the probability that the bead-spring chain has an instantaneous configuration in the range dri, . . . , drjy 
about ri, . . . , rjv, and the beads in the chain have velocities in the range rfri, . . . , drjv about f i, . . . , tn- 
The configurational distribution function can be found by integrating V over all the bead velocities, 



9(r u ...,r N ,t) = J... J VdH-..dr N (36) 
The distribution of internal configurations ip, is given by, 

i>{Q 1 ,...,Q N _ 1 ,t) = J &{r c ,Q 1 ,...,Q N _ 1 ,t)dr c (37) 

where, = ^, as a result of the Jacobian relation for the configurational vectors [ |^], 

d(ri,..., r N ) 

d{r c ,Qn Q N i ) 



= 1 



Note that the normalisation condition J tpdQ 1 dQ 2 ■ ■ ■ dQ N _ x — 1 is satisfied by ip. When the config- 
urations of the bead-spring chain do not depend on the location of the center of mass, as in the case of 
homogeneous flows with no concentration gradients, (1/V) t/j — \P, where V is the volume of the solution. 
The velocity-space distribution function 5 is defined by, 

E(n,...,r N , r 1: ...,r N , t) = — (38) 

Note that S satisfies the normalisation condition J. . .J Edri . . .dr^ — 1. Under certain circumstances 
that are discussed later, it is common to assume that the velocity-space distribution function is Maxwellian 
about the mass-average solution velocity, 



3 = Mm exp 



1 



2k B T 



m(fi — v) 2 + . . . + m(rN — v) z 



(39) 



where Mm is the normalisation constant for the Maxwellian distribution. Making this assumption implies 
that one expects the time scales involved in equilibration processes in momentum space to be much smaller 
than the time scales governing relaxation processes in configuration space. 
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Averages of quantities which are functions of the bead positions and bead velocities are defined analo- 
gously to the those in the previous section, namely, 



(X) = J...J XVdr 1 ...dr N dr 1 ...dr N (40) 
is the the phase space average of X, while the velocity- space average is, 

i X l = If XEdrx...dr N (41) 

For quantities X that depend only on the internal configurations of the polymer chain and not on the 
center of mass coordinates or bead velocities, 

(X) = J X^dQ 1 dQ 2 ...dQ N _ 1 (42) 

3.2.2. The equation of motion 

The equation of motion for a bead in a bead-spring chain is derived by considering the forces acting 
on it. The total force F^, on bead /i is, F M = ^ F^, where the F$ , i = 1,2,..., are the various 
intramolecular and solvent forces acting on the bead. The fundamental difference among the various 
molecular theories developed so far for the description of dilute polymer solutions lies in the kinds of 
forces F^' that are assumed to be acting on the beads of the chain. In almost all these theories, the 
accelaration of the beads due to the force F M is neglected. A bead-spring chain model incorporating 
bead inertia has shown that the neglect of bead inertia is justified in most practical situations [ J37|. The 
equation of motion is consequently obtained by setting F ^ = o. Here, we consider the following force 
balance on each bead fj,, 

fW+F®+F<P+fW=o (/i = l,2,...,iV) (43) 

where, F^f 1 is the hydrodynamic drag force, F^ is the Brownian force, Fjf' is the intramolecular force 

due to the potential energy of the chain, and -F„ is the force due to the presence of internal viscosity. 
These are the various forces that have been considered so far in the literature, which are believed to play 
a crucial role in determining the polymer solution's transport properties. The nature of each of these 
forces is discussed in greater detail below. Note that, as is common in most theories, external forces 
acting on the bead have been neglected. However, their inclusion is reasonably straight forward [ ^ . 

The hydrodynamic drag force F^/ , is the force of resistance offered by the solvent to the motion of the 
bead ijl. It is assumed to be proportional to the difference between the velocity-averaged bead velocity 
[[r M ]] and the local velocity of the solution, 

^ ) = -C[P>I-K + <)] ( 44 ) 

where Q is bead friction coefficient. Note that for spherical beads with radius o, in a solvent with viscosity 
r] s , the bead friction coefficient £ is given by the Stokes expression: £ = Q-Krj s a. The velocity-average of 
the bead velocity is not carried out with the Maxwellian distribution since this is just the mass-average 
solution velocity. However, it turns out that an explicit evaluation of the velocity- average is unnecessary 
for the development of the theory. Note that the velocity of the solution at bead fj, has two components, 
the imposed flow field = Vo + K(t)-r^, and the perturbation of the flow field v'^ due to the motion of the 
other beads of the chain. This perturbation is called hydrodynamic interaction, and its incorporation in 
molecular theories has proved to be of utmost importance in the prediction of transport properties. The 
presence of hydrodynamic interaction couples the motion of one bead in the chain to all the other beads, 
regardless of the distance between the beads along the length of the chain. In this sense, hydrodynamic 
interaction is a long range phenomena. 
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The perturbation to the flow field v'(r) at a point r due to the presence of a point force F(r') at the 
point r', can be found by solving the linearised Navier-Stokes equation [ fij |), 

v'(r) = fl(r - r') ■ F(r') (45) 

where O(r), called the Oseen-Burgers tensor, is the Green's function of the linearised Navier-Stokes 
equation, 

v 1 ( rr\ , i 

^=8^( 1 + ^J ™ 

The effect of hydrodynamic interaction is taken into account in polymer kinetic theory by treating the 
beads in the bead-spring chain as point particles. As a result, in response to the hydrodynamic drag 
force acting on each bead, each bead exerts an equal and opposite force on the solvent at the point that 
defines its location. The disturbance to the velocity at the bead v is the sum of the disturbances caused 
by all the other beads in the chain, v' u = — ^ fl uil {r v — r^) ■ F^\ where, — fl ufi is given by, 



1 



for /i^i/ 



n^ = { StttmvV i% ) " " (47) 

for [i = v 

The Brownian force F^\ on a bead /i is the result of the irregular collisions between the solvent 
molecules and the bead. Instead of representing the Brownian force by a randomly varying force, it is 
common in polymer kinetic theory to use an averaged Brownian force, 

F^ = -k,T(^) (48) 



As mentioned earlier, the origin of this expression can be understood within the framework of the complete 
phase space theory [ ||, || . Note that the Maxwellian distribution has been used to derive equation (Q) . 

The total potential energy <j> of the bead-spring chain is the sum of the potential energy S of the elastic 
springs, and the potential energy E due to the presence of excluded volume interactions between the 
beads. The force F^ on a bead [i due to the intramolecular potential energy <j) is given by, 

" dr 1 ' 



In addition to the various forces discussed above, the internal viscosity force F^ v \ has received consid- 
erable attention in literature [ |3|, |38|, ^] , though it appears not to have widespread acceptance. Various 
physical reasons have been cited as being responsible for the internal viscosity force. For instance, the 
hindrance to internal rotations due to the presence of energy barriers, the friction between two monomers 
on a chain that are close together in space and have a non-zero relative velocity, and so on. The simplest 
models for the internal viscosity force assume that it acts only between neighbouring beads in a bead- 
spring chain, and depends on the average relative velocities of these beads. Thus, for a bead [i that is 
not at the chain ends, 

Ff> = v ( ^-r^-r^ . f . ^ _ . ^ _ / (r M - r^) (r r^U f . ^ 
V | ^V+i -tv | / V I r M - r M _! | / 

where tp is the internal viscosity coefficient. A scaling theory for a more general model that accounts for 
internal friction between arbitrary pairs of monomers has also been developed [ |35| . 
The equation of motion for bead v can consequently be written as, 



-c 



IKJ -v -K-r u +J2^- F<W ] - k B T + F£» + F^) = ( 5 1) 
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Since = k B T (d lnf/drj - F { f> - F$ v \ equation (fl]) can be rearranged to give, 

[fJ=to + «-r» + |Eli«-( - k * T ^7 + + ) ( 52 ) 
where ■y fJiV is the dimensionless diffusion tensor [ ||, 

Ipu = <V 1 + C ^ (53) 



By manipulating equation (g2|) , it is possible to rewrite the equation of motion in terms of the velocities 
of the center of mass r c and the bead-connector vectors Q k , 



where, B kv is defined by, B^ v = b~k+\,v — $kv, the internal viscosity force, f k w \ in the direction of the 
connector vector Q k is, 

(56) 

I Qfc I 

and the tensor Aj k which accounts for the presence of hydrodynamic interaction is defined by, 

A? fc = ~ Af'fc 1 + C(^j,k + fij+l.k+l — flj,k+l — fij+l,k) (57) 

Here, Ajk is the Rouse matrix, 
r 2 for |j - fc| = 0, 

^fc = I -1 for |i-fc| = 1, (58) 
v otherwise 

In order to obtain the diffusion equation for a dilute solution of bead-spring chains, the equation of 
motion derived here must be combined with an equation of continuity. 

3.2.3. The diffusion equation 

The equation of continuity or 'probability conservation', which states that a bead-spring chain that 
disappears from one configuration must appear in another, has the form [ 

V 

The independence of \& from the location of the center of mass for homogeneous flows, and the result 
tr k = 0, for an incompressible fluid, can be shown to imply that the equation of continuity can be written 
in terms of internal coordinates alone as [ 0, 

j J 

The general diffusion equation which governs the time evolution of the instantaneous configurational 
distribution function ip, in the presence of hydrodynamic interaction, arbitrary spring and excluded 
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volume forces, and an internal viscosity force given by equation (pq), is obtained by substituting the 
equation of motion for from equation (55) into equation (60). It has the form, 



Equations such as ( |6l| ) are also referred to as Fokker-Planck or Smoluchowski equations in the literature. 
The diffusion equation ( |6l|) is the most fundamental equation of the kinetic theory of dilute polymer 
solutions since a knowledge of ip, for a flow field specified by k, would make it possible to evaluate 
averages of various configuration dependent quantities and thereby permit comparison of theoretical 
predictions with experimental observations. 

The diffusion equation can be used to derive the time evolution equation of the average of any arbitrary 
configuration dependent quantity, X(Q 1 , . . . . Qjy-i by multiplying the left and right hand sides of 
equation ( |6l| ) by X and integrating both sides over all possible configurations, 

d ( X ) ST / n dX \ kBT X^/~A dln ^ dX \ 1 Sr/A \ ^ dX \(ro\ 



dQi 



Except for a situation where nearly all the important microscopic phenomena arc neglected, the diffu- 
sion equation (|6l| ) is unfortunately in general analytically insoluble. There have been very few attempts to 
directly solve diffusion equations with the help of a numerical solution procedure [ ^, |l^| . In this context 
it is worth bearing in mind that what are usually required are averages of configuration dependent quanti- 
ties. However, in general even averages cannot be obtained exactly by solving equation (p2|) . As a result, 
it is common in most molecular theories to obtain the averages by means of various approximations. 

In order to examine the validity of these approximations it is vitally important to compare the ap- 
proximate predictions of transport properties with the exact predictions of the models. One of the ways 
by which exact numerical results may be obtained is by adopting a numerical procedure based on the 
mathematical equivalence of diffusion equations in polymer configuration space and stochastic differential 
equations for the polymer configuration [ . Instead of numerically solving the analytically intractable 
difusion equation for the distribution function, stochastic trajectories can be generated by Brownian 
dynamics simulations based on a numerical integration of the appropriate stochastic differential equa- 
tion. Averages calculated from stochastic trajectories (obtained as a solution of the stochastic differential 
equations) , are identical to the averages calculated from distribution functions (obtained as a solution of 
the diffusion equations). It has now become fairly common for any new approximate molecular theory of 
a microscopic phenomenon to establish the accuracy of the results with the help of Brownian dynamics 
simulations. In this chapter, while results of such simulations are cited, details of the development of 
the appropriate stochastic differential equations are not discussed. A comprehensive introduction to the 
development of stochastic differential equations which are equivalent to given diffusion equations for the 
probability density in configuration space, can be found in the treatise by Ottinger [ ^9f . 

3.2.4- The stress tensor 

The expression for the stress tensor in a dilute polymer solution was originally obtained by the use 
of simple physical arguments which considered the various mechanisms that contributed to the flux of 
momentum across an oriented surface in the fluid [ ^]. The major mechanisms considered were the 
transport of momentum by beads crossing the surface, and the tension in the springs that straddle the 
surface. These physical arguments help to provide an intuitive understanding of the origin of the different 
terms in the stress tensor expression. On the other hand, such arguments are difficult to pursue in the 
presence of complicated microscopic phenomena, and there is uncertainity about the completeness of the 
final expression. An alternative approach to the derivation of the expression for the stress tensor has 
been to use more fundamental arguments that consider the complete phase space of the polymeric fluid [ 

Els- 



ie 



A very general expression for the polymer contribution to the stress tensor, derived by adopting the 
complete phase space approach, for models without constraints such as the bead-spring chain model, in 
the presence of hydrodynamic interaction and an arbitrary intramolecular potential force, is the modified 
Kramers expression [ || , 

tP = u p X ( ir » - r ^ p{ - ] ) + ( N ~ 1 ) n P k * T 1 ( 63 ) 



When rewritten in terms of the internal coordinates of a bead-spring chain, equation (|63|) assumes a form 
called the Kramers expression, 

rP = -n P X<^^T> + (^- 1 )^ fe BT 1 (64) 

It is important to note that the presence of internal viscosity has not been taken into account in the 
phase space theories used to derive the modified Kramers expression (|6^). When examined from the 
standpoint of thermodynamic considerations, the proper form of the stress tensor in the presence of 
internal viscosity appears to be the Giesekus expression rather than the Kramers expression [ [39) . Since 
predictions of models with internal viscosity are not considered in this chapter, the Giesekus expression 
is not discussed here. 

In order to evaluate the stress tensor, for various choices of the potential energy <f), it turns out that it is 
usually necessary to evaluate the second moments of the bead connector vectors, ( QjQk )• An equation 
that governs the time evolution of the second moments can be obtained with the help of equation ([32]). 
It has the form, 

j t {Q 3 Qk) = K-(Q j Q k ) + {Q j Q k )- K ) + ^^(A jk ) 

- jE{<Oi[^ + /^]-^> + <^-[^ + /^]-0*>} (65) 

The second moment equation (pq), which is an ordinary differential equation, is in general not a closed 
equation for ( QjQ k ), since it invoves higher order moments on the right hand side. 

Within the context of the molecular theory developed thus far, it is clear that the prediction of the 
rheological properties of dilute polymer solutions with a bead-spring chain model usually requires the 
solution of the second moment equation ([35|). To date however, there are no solutions to the general second 
moment equation (|6^) which simultaneously incorporates the microscopic phenomena of hydrodynamic 
interaction, excluded volume, non-linear spring forces and internal viscosity. Attempts have so far been 
restricted to treating a smaller set of combinations of these phenomenon. The simplest molecular theory, 
based on a bead-spring chain model, for the prediction of the transport properties of dilute polymer 
solutions is the Rouse model. The Rouse model neglects all the microscopic phenomenon listed above, 
and consequently fails to predict many of the observed features of dilute solution behavior. In a certain 
sense, however, it provides the framework and motivation for all further improvements in the molecular 
theory. The Rouse model and its predictions are introduced below, while improvements in the treatment 
of hydrodynamic interactions alone are discussed subsequently. 

3.3. The Rouse model 

The Rouse model assumes that the springs of the bead-spring chain are governed by a Hookean spring 
force law. The only solvent-polymer interactions treated are that of hydrodynamic drag and Brownian 
bombardment. The diffusion equation ( |6l| ) with the effects of hydrodynamic interaction, excluded volume 
and internal viscosity neglected, and with a Hookean spring force law, has the form, 

dip d ( „ H v-^ „ n \ , k B T „ d dip 

-i = - E w ■ (" ■ Q* - j E ^ Q* ) + + -f- E a* Wj • ot ( fi6 ) 
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The diffusion equation (|66|) has an analytical solution since it is linear in the bead-connector vectors. 
It is satisfied by a Gaussian distribution. 

r/> (Q l5 . . . , Q N _ ± ) = M(t) exp [ - i Qi ■ (Oik ■ Q k ] (67) 

j, k 

where Af(i) is the normalisation factor, and the tensor ojfc which uniquely characterises the Gaussian 
distribution is identical to the second moment, 

a jk = (QjQk) (68) 

Note that the tensors 05 & are not symmetric, but satisfy the relation 0$% — ojy. (Further information 
on linear diffusion equations and Gaussian distributions can be obtained in the extended discussion in 
Appendix A of [ §§). 

Since the intramolecular potential in the Rouse model is only due to the presence of Hookean springs, 
it is straight forward to see that the Kramers expression for the stress tensor r p , is given by 

t p = -n p H °3j + (N-1) n p k B T 1 (69) 
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The tensors 05 j are obtained by solving the second moment equation (|6^), which becomes a closed 
equation for the second moments when the Rouse assumptions are made. It has the form, 

d j- 2/cbT H v— * r -1 r*7c\\ 

— (Tjk — K ■ Ujk — CTj k ■ K = — - — Ajk 1 — [ ^jm A m k + Aj m (T mk \ (70) 



Note that the solution of equation (|70|) also leads to the complete specification of the Gaussian configu- 
rational distribution function ip. 

A Hookean dumbbell model, which is the simplest example of a bead-spring chain, is obtained by 
setting N = 2. It is often used for preliminary calculations since its simplicity makes it possible to obtain 
analytical solutions where numerical solutions are unavoidable for longer chains. For such a model, 
substituting for a\i in terms of r p from equation (^) into equation (|70|), leads to following equation for 
the polymer contribution to the stress tensor, 

r p + X H t p {1) = -n p k B TX H 7 (71) 

where = dr p jdt — k ■ r p — t p ■ , is the convected time derivative [ [| of t p , and Xh = (C,/AH) is 
a time constant. Equation ( [7l| ) indicates that a Hookean dumbbell model with the Rouse assumptions 
leads to a convected Jeffreys model or Oldroyd-B model as the constitutive equation for a dilute polymer 
solution. This is perhaps the simplest coarse-grained microscopic model capable of reproducing some of 
the macroscopic rheological properties of dilute polymer solutions. 

In the case of a bead-spring chain with N > 2, it is possible to obtain a similar insight into the nature 
of the stress tensor by introducing normal coordinates. These coordinates help to decouple the connector 
vectors Q 1 , . . . , Q N _i, which are coupled to each other because of the Rouse matrix. 

The connector vectors are mapped to a new set of normal coordinates, . . . , Qjv_i with the trans- 
formation, 

Qj = E U ^Q' k (72) 

k 

where, Ujk are the elements of an orthogonal matrix with the property 

(rr 1 )^ = rife,, such that, E n mj -n m/c = S jk (73) 
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The orthogonal matrix Hjk, which will henceforth be referred to as the Rouse orthogonal matrix, 
diagonalises the Rouse matrix Ajk, 

UjiAjklLki = aiSu (74) 

j, k 

where, the Rouse eigenvalues ai are given by a; = 4 sin 2 (In /2N). The elements of the Rouse orthogonal 
matrix are given by the expression, 

The diffusion equation in terms of these normal coordinates, admits a solution for the configurational 
distribution function of the form [ ^ 

AT-l 

^(Qi,...,Q / JV _ 1 )= n ^ 
fe=i 

As a consequence, the diffusion equation becomes uncoupled and can be simplified to (N — 1) diffusion 
equations, one for each of the ipk (Qk)- Since the Q' k are independent of each other, all the covariances 
(Q'jQ'k) with j ^ k are zero, and only the (N — 1) variances = (QjQj) are non-zero. Evolution 
equations for the variances o^- can then be derived from these uncoupled diffusion equations with the help 
of a procedure similar to that used for the derivation of equation ( |65| ) . 
The stress tensor is given in terms of by the expression, 

- p = E ( 77 ) 

where, 

t v 3 = -n p H cfj + n p k B T l (78) 



On substituting for o^- in terms of from equation flTq ) into the evolution equation for o^-, one obtains, 

T* j +* j T P H1) = -n P k B T\ j <y (79) 

where, the relaxation times are given by Xj = (£/2_ff a,j). Consequently, each of the satisfy 
an equation identical to equation (|7l|) for the polymer contribution to the stress tensor in a Hookean 
dumbbell model. The Rouse model, therefore, leads to a constitutive equation that is a multimode 
generalization of the convected Jeffreys or Oldroyd B model. 

ft is clear from above discussion that the process of transforming to normal coordinates enables one to 
derive a closed form expression for the stress tensor, and to gain the insight that the Rouse chain with 
beads has N independent relaxation times which describe the different relaxation processes in the chain, 
from the entire chain to successively smaller sub-chains, ft is straight forward to show that for large N, 
the longest relaxation times Xj scale with chain length as N 2 . 

A few important transport property predictions which show the limitations of the Rouse model are 
considered briefly below, ft is worth noting that since the Rouse model does not include the effect of 
excluded volume, its predictions are restricted to dilute solutions of polymers in theta solvents. This 
restriction is infact applicable to all the models of hydrodynamic interaction treated here. 

In steady simple shear flow, with n(t) given by equation (^), the three independent material functions 
that characterise such flows are f , 



MT^Xj] *i = 2n p k B Tj2rf; *2 = (80) 
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It is clear that the Rouse model accounts for the presence of viscoelasticity through the prediction of a 
nonzero first normal stress difference in simple shear flow. However, it does not predict the nonvanishing 
of the second normal stress difference, and the shear rate dependence of the viscometric functions. 

From the definition of intrinsic viscosity ([!]) and the fact that p p ~ N n p , it follows from equation ( Mj ) 
that for the Rouse model, [r)]o ~ N. This is at variance with the experimental results discussed earlier, 
and displayed in equation It is also straight forward to see that the Rouse model predicts that the 
characteristic relaxation time scales as the square of the chain length, X p ~ iV 2 . 

In small amplitude oscillatory shear, n(t) is given by equation (ph, and expressions for the material 
functions G' and G" in terms of the relaxation times Xj can be easily obtained [ || . In the intermediate 
frequency range, where as discussed earlier, experimental results indicate that both G' and G" scale as 
w 2 / 3 , the Rouse model predicts a scaling a; 1 / 2 [0. 

The translational diffusion coefficient D for a bead-spring chain at equilibrium can be obtained by 
finding the average friction coefficient Z for the entire chain in a quiescent solution, and subsequently 
using the Nernst-Einstein equation, D = k^T Z^ 1 [ p). It can be shown that for the Rouse model 
Z = C N, ie. the total friction coefficient of the chain is a sum of the individual bead friction coefficients. 
As a result, the Rouse model predicts that the diffusion coefficient scales as the inverse of the molecular 
weight. This is not observed in dilute solutions. Instead experiments indicate the scaling depicted in 
equation (|l2]). 

The serious shortcomings of the Rouse model highlighted above have been the motivation for the 
development of more refined molecular theories. The scope of this chapter is restricted to reviewing 
recent advances in the treatment of hydrodynamic interaction. 



4. HYDRODYNAMIC INTERACTION 

Hydrodynamic interaction, as pointed out earlier, is a long range interaction between the beads which 
arises because of the solvent's capacity to propagate one bead's motion to another through perturbations 
in its velocity field. It was first introduced into framework of polymer kinetic theory by Kirkwood 
and Riseman [ [l4|]. As we have seen in the development of the general diffusion equation above, it is 
reasonably straight forward to include hydrodynamic interaction into the framework of the molecular 
theory. However, it renders the resultant equations analytically intractable and as a result, various 
attempts have been made in order to solve them approximately. 

In this section, we review the various approximation schemes introduced over the years. The primary 
test of an approximate model is of course its capacity to predict experimental observations. The accuracy 
of the approximation however, can only be assessed by checking the proximity of the approximate results 
to the exact numerical results obtained by Brownian dynamics simulations. Finally, the usefulness of an 
approximation depends on its computational intensity. The individual features and deficiencies of the 
different approximations will be examined in the light of these observations. 

In the presence of hydrodynamic interaction, and with excluded volume and internal viscosity neglected, 
a bead-spring chain with Hookean springs has a configurational distribution function tp that must satisfy 
the following simplified form of the diffusion equation ([H]), 
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while the second moment equation (|65|) assumes the form, 



j t (QjQk) = « • (QjQk) + (QjQk) ■ ^ + ^f- (A jh ) - f E 
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(82) 



Equation ( p2j ) is not a closed equation for the second moments since it involves more complicated mo- 
ments on the right hand side. This is the central problem of all molecular theories which attempt to 
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predict the rheological properties of dilute polymer solutions and that incorporate hydrodynamic inter- 
action. The different approximate treatments of hydrodynamic interaction, which are discussed roughly 
chronologically below, basically reduce to finding a suitable closure approximation for the second moment 
equation. 

4.1. The Zimm model 

The Zimm model was the first attempt at improving the Rouse model by introducing the effect of hy- 
drodynamic interaction in a preaveraged or equilibrium-averaged form. The preaveraging approximation 
has been very frequently used in polymer literature since its introduction by Kirkwood and Riseman [ [l4| . 
The approximation consists of evaluating the average of the hydrodynamic tensor with the equilibrium 
distribution function (|34|), and replacing the hydrodynamic interaction tensor Ajk, wherever it occurs in 
the governing equations, with its equilibrium average Ajk- (Note that the incorporation of the effect of 
hydrodynamic interaction does not alter the equilibrium distribution function, which is still given by ( |34| ) 
for bead-spring chains with Hookean springs.) The matrix Ajk is called the modified Rouse matrix, and 
is given by, 

A ik = A lk + V2h*( ■ 2 1 1 | (83) 



where, h* = aW (H/nksT) is the hydrodynamic interaction parameter. The hydrodynamic interaction 
parameter is approximately equal to the ratio of the bead radius to the equilibrium root mean square 
length of a single spring of the bead-spring chain. This implies that h* < 0.5, since the beads cannot 
overlap. Typical values used for h* are in the range 0.1 < h* < 0.3 [ g2| . 

By including the hydrodynamic interaction in an averaged form, the diffusion equation remains linear 
in the connector vectors, and consequently is satisfied by a Gaussian distribution ( |67| ) as in the Rouse 
case. However, the covariance tensors <Tjk are now governed by the set of differential equations ( |70| ) with 
the Rouse matrix Ajk replaced with the modified Rouse matrix Ajk- Note that this modified second 
moment equation is also a closed set of equations for the second moments. 

As in the Rouse case, it is possible to simplify the solution of the Zimm model by carrying out a 
diagonalisation procedure. This is achieved by mapping the connector vectors to normal coordinates, as 
in (|72|), but in this case the Zimm orthogonal matrix H^, which diagonalises the modified Rouse matrix, 

ILjiAjkUki = aiSa (84) 

j, k 

must be found numerically for N > 4. Here, ai are the so called Zimm eigenvalues. The result of this 
procedure is to render the diffusion equation solvable by the method of separation of variables. Thus, as 
in the Rouse case, only the (N — 1) transformed coordinate variances dj are non-zero, and differential 
equations governing these variances can be derived by manipulating the uncoupled diffusion equations. 

The diagonalisation procedure enables the polymer contribution to the stress tensor t p in the Zimm 
model to be expressed as a sum of partial stresses as in equation (|77]), but the t v ^ now satisfy 
equation (|79| ) with the 'Rouse' relaxation times Xj replaced with 'Zimm' relaxation times Xj. The Zimm 
relaxation times are defined by Xj = (£/2Ha,j). 

From the discussion above, it is clear that the Zimm model differs from the Rouse model only in the 
spectrum of relaxation times. As we shall see shortly, this leads to a significant improvement in the 
prediction of linear viscoelastic properties and the scaling of transport properties with molecular weight 
in theta solvents. The Zimm model therefore establishes unequivocally the importance of the microscopic 
phenomenon of hydrodynamic interaction. On the other hand, it does not lead to any improvement in the 
prediction of nonlinear properties, and consequently subsequent treatments of hydrodynamic interaction 
have concentrated on improving this aspect of the Zimm model. 
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By considering the long chain limit of the Zimm model, ie., A — > oo, it is possible to discuss the 
universal properties predicted by the model. The various power law dependences of transport properties 
on molecular weight, characterised by universal exponents, and universal ratios formed from the prefactors 
of these dependences can be obtained. These predictions are ideal for comparison with experimental data 
on high molecular weight polymer solutions since they are parameter free. We shall discuss some universal 
exponents predicted by the Zimm model below, while universal ratios are discussed later in the chapter. 

As mentioned above, the first noticeable change upon the introduction of hydrodynamic interaction is 
the change in the relaxation spectrum. In the long chain limit, the longest relaxation times Xj scale with 
chain length as A 3 / 2 [ ^9), whereas we had found earlier that the chain length dependence of the longest 
relaxation times in the Rouse model was A 2 . 

In steady simple shear flow, the Zimm model like the Rouse model, fails to predict the experimentally 
observed occurance of non-zero second normal stress differences and the experimentally observed shear 
rate dependence of the viscometric functions. It does however lead to an improved prediction of the 
scaling of the zero shear rate intrinsic viscosity with molecular weight, [t)]q ~ A 1 / 2 . This prediction is in 
agreement with experimental results for the Mark-Houwink exponent in theta solvents (see equation (|Tl|)). 
As with the longest relaxation times, the characteristic relaxation time A p ~ A 3 / 2 . 

In small amplitude oscillatory shear, the Zimm model predicts that the material functions G and G" 
scale with frequency as u 2 / 3 in the intermediate frequency range. This is in exceedingly good agreement 
with experimental results [ |[ [lTj . 

The translational diffusion coefficient D for chainlike molecules at equilibrium, with preaveraged hydro- 
dynamic interaction, was originally obtained by Kirkwood [ [lij . Subsequently, several workers obtained a 
correction to the Kirkwood diffusion coefficient for the Zimm model [ . The exact results differ by less 
than 2% from the Kirkwood value for all values of the chain length and h* . Interestingly, three different 
approaches to obtaining the diffusion coefficient, namely, the Nernst-Einstein equation, the calculation 
of the mean-square displacement caused by Brownian forces, and the study of the time evolution of con- 



centration gradients, lead to identical expressions for the diffusion coefficient [ 21 . In the limit of very 
long chains, it can be shown that D ~ A" 1 / 2 . The Zimm model therefore gives the correct dependence 
of translational diffusivity on molecular weight in theta solvents. 

The Zimm result for the translational diffusivity has been traditionally interpreted to mean that the 
polymer coil in a theta solvent behaves like a rigid sphere, with radius equal to the root mean square 
end-to-end distance. This follows from the fact that the diffusion coefficient for a rigid sphere scales as 
the inverse of the radius of the sphere, and in a theta solvent, ( r 2 ^ scales with chain length as A. The 
solvent inside the coil is believed to be dragged along with the coil, and the inner most beads of the bead- 
spring chain are considered to be shielded from the velocity field due to the presence of hydrodynamic 
interaction [ 44, [l?]]. This intuitive notion has been used to point out the difference between the Zimm 
and the Rouse model, where all the A beads of the polymer chain arc considered to be exposed to the 
applied velocity field. Recently, by explicitly calculating the velocity field inside a polymer coil in the 
Zimm model, Ottinger [ |3^] has shown that the solvent motion inside a polymer coil is different from 
that of a rigid sphere throughout the polymer coil, and that shielding from the velocity field occurs only 
to a certain extent. 



4.2. The consistent averaging approximation 

The first predictions of shear thinning were obtained when hydrodynamic interaction was treated in a 
more precise manner than that of preaveraging the hydrodynamic interaction tensor. In order to make the 
diffusion equation ( |Sl| ) linear in the connector vectors, as pointed out earlier, it is necessary to average the 
hydrodynamic interaction tensor. However, it is not necessary to preaverage the hydrodynamic interaction 
tensor with the equilibrium distribution. On the other hand, the average can be carried out with the non- 
equilibrium distribution function ( p7|). The linearity of the diffusion equation ensures that its solution is a 
Gaussian distribution. Ottinger [ |20|, |22] | suggested that the hydrodynamic interaction tensor occuring 
in the diffusion equation be replaced with its non-equilibrium average. Since it is necessary to know 
the averaged hydrodynamic interaction tensor in order to find the non-equilibrium distribution function, 
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both the averaged hydro-dynamic interaction tensor and the non-equilibrium distribution function must 
be obtained in a self- consistent manner. 

Several years ago, Fixman [ ^| introduced an iterative scheme (beginning with the equilibrium dis- 
tribution function), for refining the distribution function with which to carry out the average of the 
hydrodynamic interaction tensor. The self-consistent scheme of Ottinger is recovered if the iterative 
procedure is repeated an infinite number of times. However, Fixman carried out the iteration only upto 
one order higher than the preaveraging stage. 

The average of the hydrodynamic interaction tensor evaluated with the Gaussian distribution ( |67| ) is 
an (N — 1) x (N — 1) matrix with tensor components, Ajk, defined by, 



Aj k = A jk 1 + V2h 
where the tensors are given by, 
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and the function of the second moments, H(cr) is, 

*w = ^ / dk h ( 1 - § ) • fe) (87) 

Note that the convention H(djj)/0 — has been adopted in equation ( |85| ) above. 

The self-consistent closure approximation therefore consists of replacing the hydrodynamic interaction 
tensor Ajk in equation ( p2[ ) with its non-equilibrium average Ajk- As in the earlier approximations, 
this leads to a system of (N — l) 2 coupled ordinary differential equations for the components of the 
covariance matrix ajk- Their solution permits the evaluation of the stress tensor through the Kramers 
expression ([39|), and as a consequence all the relevant material functions. 

Viscometric functions in steady simple shear flows were obtained by Ottinger [ for chains with 
N < 25 beads, while material functions in start-up of steady shear flow, cessation of steady shear flow, 
and stress relaxation after step-strain were obtained by Wedgewood and Ottinger [ for chains with 
N < 15 beads. The latter authors also include consistently-averaged FENE springs in their model. 

Shear rate dependent viscometric functions, and a nonzero positive second normal stress difference are 
predicted by the self-consistent averaging approximation; a marked improvement over the predictions 
of the Zimm model. Both the reduced viscosity and the reduced first normal stress difference initially 
decrease with increasing shear rate. However, for long enough chains, they begin to rise monotonically at 
higher values of the reduced shear rate j3. This rise is a consequence of the weakening of hydrodynamic 
interaction in strong flows due to an increase in the separation between the beads of the chain. With 
increasing shear rate, the material functions tend to the shear rate independent Rouse values, which (for 
long enough chains), are higher than the zero shear rate consistently-averaged values. The prediction 
of shear thickening behavior is not in agreement with the shear thinning that is commonly observed 
experimentally. However, as mentioned earlier, some experiments with very high molecular weight systems 
seem to suggest the existence of shear thinning followed by shear thickening followed again by shear 
thinning as the shear rate is increased. While only shear thickening at high shear rates is predicted with 
Hookean springs, the inclusion of consistently-averaged FENE springs in the model leads to predictions 
which are qualitatively in agreement with these observations, with the FENE force becoming responsible 
for the shear thinning at very high shear rates [ |d], |l5| . 

The means of examining the accuracy of various approximate treatments of hydrodynamic interaction 
was established when the problem was solved exactly with the help of Brownian dynamics simulations 
with full hydrodynamic interaction included [ H6], Bq] . These simulations reveal that while the predictions 
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of the shear rate dependence of the viscosity and first normal stress difference by the self-consistent 
averaging procedure are in qualitative agreement with the Brownian dynamics simulations, they do not 
agree quantitatively. Further, in contrast to the consistent-averaging prediction, at low shear rates, a 
negative value for the second normal stress difference is obtained. As noted earlier, the sign of the second 
normal stress difference has not been conclusively established [ [| . 

The computational intensity of the consistent-averaging approximation leads to an upper bound on 
the length of chain that can be examined. As a result, it is not possible to discuss the universal shear 
rate dependence of the viscometric functions predicted by it. On the other hand, it is possible to come to 
certain general conclusions regarding the nature of the stress tensor in the long chain limit, and to predict 
the zero shear rate limit of certain universal ratios [pp|j. Thus, it is possible to show the important result 
that the polymer contribution to the stress tensor depends only on a length scale and a time scale, and 
not on the strength of the hydrodynamic interaction parameter h*. In the long chain limit, h* can be 
absorbed into the basic time constant, and it does not occur in any of the non-dimensional ratios. Indeed 
this is also true of the finite extensibility parameter 6, which can also be shown to have no influence on 
the long chain rheological properties [ [22) . The long chain limit of the consistent-averaging approximation 
is therefore a parameter free model. 

It is possible to obtain an explicit representation of the modified Kramers matrix for infinitely long 
chains by introducing continuous variables in place of discrete indices [ . This enables the analytical 
calculation of various universal ratios predicted by the consistent-averaging approximation. These predic- 
tions are discussed later in this chapter. However, two results are worth highlighting here. Firstly, it can 
be shown explicitly that the leading order corrections to the large N limit of the various universal ratios 
are of order (1/y/N), and secondly, there is a special value of h* — 0.2424..., at which the leading order 
corrections are of order (1/N). These results have proven to be very useful for subsequent numerical 
exploration of the long chain limit in more accurate models of the hydrodynamic interaction. 

Short chains with consistently-averaged hydrodynamic interaction, as noted earlier, do not show shear 
thickening behavior; this aspect is revealed only with increasing chain length. Furthermore, it is not clear 
with the kind of chain lengths that can be examined, whether the minimum in the viscosity and first 
normal stress curves continue to exist in the long chain limit [ ^| . The examination of long chain behavior 
is therefore important since aspects of polymer solution behavior might be revealed that are otherwise 
hidden when only short chains are considered. The introduction of the decoupling approximation by 
Magda, Larson and Mackay [ |T^] and Kishbaugh and McHugh [ |l5| made the examination of the shear 
rate dependence of long chains feasible. The decoupling approximation retains the accuracy of the self- 
consistent averaging procedure, but is much more computationally efficient. 



4.3. The decoupling approximation 

The decoupling approximation introduced by Magda et al. [ [ll| and Kishbaugh and McHugh [ (who 
use FENE springs in place of the Hookean springs of Magda et al.) consists of extending the 'diagonalisc 
and decouple' procedure of the Rouse and Zimm theories to the case of the self consistently averaged 
theory. They first transform the connector vectors Qj to a new set of coordinates Q'j using the time- 
invariant Rouse orthogonal matrix Hjk- (Kishbaugh and McHugh also use the Zimm orthogonal matrix). 
The same orthogonal matrix H^ is then assumed to diagonalise the matrix of tensor components Ajk- 
While the process of diagonalisation was exact in the Rouse and Zimm theories, it is an approximation 
in the case of the decoupling approximation. It implies that even in the self consistently averaged theory 
the diffusion equation can be solved by the method of separation of variables, and only the (N — 1) 
transformed coordinate variances dj are non-zero. The differential equations governing these variances 
can be derived from the uncoupled diffusion equations and solved numerically. The appropriate material 
functions are then obtained using the Kramers expression in terms of the transformed coordinates, namely, 
equations @ and J78| ). 

The decrease in the number of differential equations to be solved, from (N — l) 2 for the covariances 
ajk to (N — 1) for the variances e^-, is suggested by Kishbaugh and McHugh as the reason for the great 



reduction in computational time achieved by the decoupling approximation. Prakash and Ottinger [ 33 
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discuss the reasons why this argument is incomplete, and point out the inconsistencies in the decoupling- 
procedure. Furthermore, since the results are only as accurate as the consistent averaging approximation, 
the decoupling approximation is not superior to the consistent averaging method. However, these papers 
are important since the means by which a reduction in computational intensity may be achieved, without 
any significant sacrifice in accuracy, was first proposed in them. Further, the persistence of the minimum 
in the viscosity and first normal stress curves even for very long chains, and the necessity of including 
FENE springs in order to generate predictions in qualitative agreement with experimental observations 
in high molecular weight systems, is clearly elucidated in these papers. 

4.4. The Gaussian approximation 

The closure problem for the second moment equation is solved in the preaveraging assumption of Zimm, 
and in the self consistent averaging method of Ottinger, by replacing the tensor Aj k with an average. 
As a result, fluctuations in the hydrodynamic interaction are neglected. The Gaussian approximation [ 
p4j Eq, E4 makes no assumption with regard to the hydrodynamic interaction, but assumes that the 



solution of the diffusion equation (81) may be approximated by a Gaussian distribution (J67|). Since all 
the complicated averages on the right hand side of the second moment equation ([32j) can be reduced to 
functions of the second moment with the help of the Gaussian distribution, this approximation makes it 
a closed equation for the second moments. 

The evolution equation for the covariances (Tjk is given by, 
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L otherwise 

The convention K(tTjj)/Q — 0, has been adopted in equation (|S9|). Both the hydrodynamic interaction 
functions H{a) and K{a) can be evaluated analytically in terms of elliptic integrals. The properties of 
these functions are discussed in great detail in the papers by Ottinger and coworkers [ ^2|, |2^, [f8|, [l?J . 

All the approximations discussed earlier (with the exception of the decoupling approximation) can be 
derived by a process of succesive simplification of the explicit results for the Gaussian approximation given 
above. The equations that govern the self consistently averaged theory can be obtained by dropping the 
last term in equation (|88|), which accounts for the presence of fluctuations. Replacing H(tr) by 1 in these 
truncated equations leads to the governing equations of the Zimm model, while setting h* = leads to 
the Rouse model. 
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Material functions predicted by the Gaussian approximation in any arbitrary homogeneous flow may 
be obtained by solving the system of (N — l) 2 coupled ordinary differential equations for the components 
of the covariance matrix <Tjk (88). Small amplitude oscillatory shear flows and steady shear flow in the 
limit of zero shear rate have been examined by Ottinger [ |24| for chains with N < 30 beads, while 
Zylka [ [l8| has obtained the material functions in steady shear flow for chains with N < 15 beads and 
compared his results with those of Brownian dynamics simulations (the comparison was made for chains 
with N = 12 beads). 

The curves predicted by the Gaussian approximation for the storage and loss modulus, G and G", 
as a function of the frequency to, are nearly indistinguishable from the predictions of the Zimm theory, 
suggesting that the Zimm approximation is quite adequate for the prediction of linear visco-elastic prop- 
erties. There is, however, a significant difference in the prediction of the relaxation spectrum. While 
the Zimm model predicts a set of (N — 1) relaxation times with equal relaxation weights, the Gaussian 
approximation predicts a much larger set of relaxation times than the number of springs in the chain, with 
relaxation weights that are different and dependent on the strength of the hydrodynamic interaction [ ^8| . 
These results indicate that entirely different relaxation spectrum lead to similar curves for G' and G", 
and calls into question the common practice of obtaining the relaxation spectrum from experimentally 
measured curves for G' and G" (see also the discussion in [ Q). 

The zero shear rate viscosity and first normal stress difference predicted by the Gaussian approximation 
are found to be smaller than the Zimm predictions for all chain lengths. By extrapolating finite chain 
length results to the infinite chain limit, Ottinger has shown that this reduction is by a factor of 72% 
- 73%, independent of the strength of the hydrodynamic interaction parameter. Other universal ratios 
predicted by the Gaussian approximation in the limit of zero shear rate are discussed later in the chapter. 

A comparison of the predicted shear rate dependence of material functions in simple shear flow with the 
results of Brownian dynamics simulations reveals that of all the approximate treatments of hydrodynamic 
interaction introduced so far, the Gaussian approximation is the most accurate [ ^|. Indeed, at low 
shear rates, the negative second normal stress difference predicted by the Gaussian approximation is in 
accordance with the simulations results. 

Inspite of the accuracy of the Gaussian approximation, its main drawback is its computational intensity, 
which renders it difficult to examine chains with large values of N. Apart from the need to examine long 
chains for the reason cited earlier, it also necessary to do so in order to obtain the universal predictions 
of the model. A recently introduced approximation which enables the evaluation of universal viscometric 
functions in shear flow is discussed in the section below. Before doing so, however, we first discuss the 
significant difference that a refined treatment of hydrodynamic interaction makes to the prediction of 
translational diffusivity in dilute polymer solutions. 

The correct prediction of the scaling of the diffusion coefficient with molecular weight upon introduction 
of pre-averaged hydrodynamic interaction in the Zimm model demonstrates the significant influence 
that hydrodynamic interaction has on the translational diffusivity of the macromolecule. While the 
pre-averaging assumption appears adequate at equilibrium, it predicts a shear rate independent scalar 
diffusivity even in the presence of a flow field. On the other hand, both the improved treatments of 
hydrodynamic interaction, namely, consistent averaging and the Gaussian approximation, reveal that 
the translational diffusivity of a Hookean dumbbell in a flowing homogeneous solution is described by 
an anisotropic diffusion tensor which is flow rate dependent [ ^3[ ^7). Indeed, unlike in the Zimm 
case, the three different approaches mentioned earlier for calculating the translational diffusivity do not 
lead to identical expressions for the diffusion tensor [ ^3[ |2q] . Insight into the origin of the anisotropic 
and flow rate dependent behavior of the translational diffusivity is obtained when the link between the 
polymer diffusivity and the shape of the polymer molecule in flow [ |l2| is explored [ |3l], |32|]. It is 
found that the solvent flow field alters the distribution of mass about the centre of the dumbbell. As a 
consequence, the dumbbell experiences an average friction that is anisotropic and flow rate dependent. 
The discussion of the influence of improved treatments of hydrodynamic interaction on the translational 
diffusivity has so far been confined to the Hookean dumbbell model. This is because the concept of the 
center of resistance, which is very useful for simplifying calculations for bead-spring chains in the Zimm 
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case, cannot be employed in these improved treatments [ ^3) . 

4.5. The twofold normal approximation 

The twofold normal approximation borrows ideas from the decoupling approximation of Magda et al. [ 
p"8[ and Kishbaugh and McHugh [ [l5| in order to reduce the computational intensity of the Gaussian 
approximation. As in the case of the Gaussian approximation, and unlike in the case of the consistent- 
averaging and decoupling approximations where it is neglected, fluctuations in the hydrodynamic inter- 
action are included. In a sense, the twofold normal approximation is to the Gaussian approximation, 
what the decoupling approximation is to the consistent-averaging approximation. The computational 
efficiency of the decoupling approximation is due both to the reduction in the set of differential equations 
that must be solved in order to obtain the stress tensor, and to the procedure that is adopted to solve 
them [ These aspects are also responsible for the computational efficiency of the twofold normal 
approximation. However, the derivation of the reduced set of equations in the twofold normal approx- 
imation is significantly different from the scheme adopted in the decoupling approximation; it is more 
straight forward, and avoids the inconsistencies that are present in the decoupling approximation. 

Essentially the twofold normal approximation, (a) assumes that the configurational distribution func- 
tion ip is Gaussian, (b) uses the Rouse or the Zimm orthogonal matrix Ujk to map Qj to 'normal' 
coordinates Q'j, and (c) assumes that the covariance matrix <Tjk is diagonalised by the same orthogonal 
matrix, ie. J2j k ^jp °5fc Hfcg = (Q' P Q' q ) = <W This leads to the following equations for the (N — 1) 
variances e^ , 
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and the matrix is given by, 
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In equations (^3|) and (|94|), the tensors Ajk and Ti p , mn are given by equations (85) and (^9|), respectively. 
However, the argument of the hydrodynamic interaction functions is now given by, 
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The decoupling approximation is recovered from the twofold normal approximation when the last term 
in equation (92), which accounts for fluctuations in hydrodynamic interaction, is dropped. Thus the two 



different routes for finding governing equations for the quantities lead to the same result. However, 
Prakash and Ottinger [ ^3) have shown that this is in some sense a fortuitous result, and indeed the key 
assumption made in the decoupling approximation regarding the diagonalisation of Ajk is not tenable. 
The Zimm model in terms of normal modes may be obtained from equation ( [92] ) by dropping the last 
term, and substituting Ajk in place of Ajk- Of course the Zimm orthogonal matrix must be used to 
carry out the diagonalisation in equation (B3) . The diagonalised Zimm model reduces to the diagonalised 
Rouse model upon using the Rouse orthogonal matrix and on setting h* = 0. 

The evolution equations (|9^) have been solved to obtain the zero shear rate properties for chains with 
N < 150, when the Zimm orthogonal matrix is used for the purpose of diagonalisation, and for chains 
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Figure 1. Universal viscometric functions in theta solvents. Reproduced from [ p3[ . 

with N < 400, when the Rouse orthogonal matrix is used. Viscometric functions at finite shear rates in 
simple shear flows have been obtained for chains with N < 100 [ |33| . The results are very close to those 
of the Gaussian approximation; this implies that they must also lie close to the results of exact Brownian 
dynamics simulations. The reasons for the reduction in computational intensity of the twofold normal 
approximation are discussed in some detail in [ |3i| . The most important consequence of introducing the 
twofold normal approximation is that rheological data accumulated for chains with as many as 100 beads 
can be extrapolated to the limit N — > oo, and as a result, universal predictions may be obtained. 

4.6. Universal properties in theta solvents 

One of the most important goals of examining the influence of hydrodynamic interactions on polymer 
dynamics in dilute solutions is the calculation of universal ratios and master curves. These properties do 
not depend on the mechanical model used to represent the polymer molecule. Consequently, they reflect 
the most general consequence of the way in which hydrodynamic interaction has been treated in the 
theory. They are also the best means to compare theoretical predictions with experimental observations 
since they are parameter free. 

There appear to be two routes by which the universal predictions of models with hydrodynamic in- 
teraction have been obtained so far, namely, by extrapolating finite chain length results to the limit of 
infinite chain length where the model predictions become parameter free, and by using renormalisation 
group theory methods. 

In the former method, there are two essential requirements. The first is that rheological data for 
finite chains must be generated for large enough values of N so as to be able to extrapolate reliably, ie. 
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Table 1 

Universal ratios in the limit of zero shear rate. The exact Zimm values and the Gaussian approximation 
(GA) values for U v r and Uy v are reproduced from []2{|, the exact consistent-averaging values from [ p0[ , 
and the renormalisation group (RG) results from [ E5| . The twofold normal approximation values with 
the Zimm orthogonal matrix (TNZ) and the remaining GA values are reproduced from [|3^]. Numbers 
in parentheses indicate the uncertainity in the last figure. 













Zimm 


2.39 


1.66425 


0.413865 


0.0 


CA 


2.39 


1.66425 


0.413865 


0.010628 


GA 


1.835 (1) 


1.213(3) 


0.560 (3) 


-0.0226 (5) 


RG 




1.377 


0.6096 


-0.0130 


TNZ 


1.835 (1) 


1.210(2) 


0.5615(3) 


-0.0232(1) 



with small enough error, to the limit N — > oo. The second is that some knowledge of the leading order 
corrections to the infinite chain length limit must be obtained in order to carry out the extrapolation in an 
efficient manner. It is clear from the discussion of the various approximate treatments of hydrodynamic 
interaction above that it is possible to obtain universal ratios in the zero shear rate limit in all the cases. 
Four universal ratios that are frequently used to represent the rheological behavior of dilute polymer 
solutions in the limit of zero shear rate are [ ^9) , 

UrlX - nfcsTAr ^ „™ nr]s (4^/3) 



— 2 = — — 

%,o *i,o 



where, Ai is the longest relaxation time, and R g is the root-mean-square radius of gyration at equilibrium. 
With regard to the leading order corrections to these ratios, it has been possible to obtain them explicitly 
only in the consistently-averaged case [ [2^] . In both the Gaussian approximation and the twofold normal 
approximation it is assumed that the leading order corrections are of the same order, and extrapolation 
is carried out numerically by plotting the data as a function of (1/y/N). Because of their computational 
intensity, it is not possible to to obtain the universal shear rate dependence of the viscometric functions 
predicted by the consistent-averaging and Gaussian approximations. However, it is possible to obtain 
these master curves with the twofold normal approximation. 

Table |l| presents the prediction of the universal ratios (|96| ) by the various approximate treatments. 
Miyaki et al. [ [l9| have experimentally obtained a value of U n u — 1.49 (6) for polystyrene in cyclohexane 
at the theta temperature. Figure [l] displays the viscometric functions predicted by the two fold normal 
approximation. The coincidence of the curves for the different values of h* indicate the parameter free 
nature of these results. Divergence of the curves at high shear rates implies that the data accumulated 
for chains with N < 100 is insufficient to carry out an accurate extrapolation at these shear rates. The 
incorporation of the effect of hydrodynamic interaction into kinetic theory clearly leads to the prediction 
of shear thickening at high shear rates even in the long chain limit. 

In both table |l| and figure [l], the results of renormalisation group calculations (RG) are also presented [ 
p^ , [i"7| . As mentioned earlier, the renormalisation group theory approach is an alternative procedure for 
obtaining universal results. It is essentially a method for refining the results of a low-order perturbative 
treatment of hydrodynamic interaction by introducing higher order effects so as to remove the ambiguous 
definition of the bead size. All the infinitely many interactions for long chains are brought in through the 
idea of self-similarity. It is a very useful procedure by which a low-order perturbation result, which can 
account for only a few interactions, is turned into something meaningful. However, systematic results can 
only be obtained near four dimensions, and one cannot estimate the errors in three dimensions reliably. 
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The Gaussian and twofold normal approximations on the other hand are non-perturbative in nature, and 
are essentially 'uncontrolled' approximations with an infinite number of higher order terms. 

It is clear from the figures that the two methods lead to significantly different results at moderate to 
high shear rates. A minimum in the viscosity and first normal stress difference curves is not predicted by 
the renormalisation group calculation, while the twofold normal approximation predicts a small decrease 
from the zero shear rate value before the monotonic increase at higher shear rates. The good comparison 
with the results of Brownian dynamics simulations for short chains indicates that the twofold normal 
approximation is likely to be more accurate than the renormalisation group calculations. 

5. CONCLUSIONS 

This chapter discusses the development of a unified basis for the treatment of non-linear microscopic 
phenomena in molecular theories of dilute polymer solutions and reviews the recent advances in the 
treatment of hydrodynamic interaction. In particular, the successive refinements which ultimately lead 
to the prediction of universal viscometric functions in theta solvents have been highlighted. 
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